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Introduction In the January 2006 issue of The College Mathematics Journal, M. 
Andreoli posed the following problem [1]: 



This article is written with the hope to draw attention to a method that has 
been proposed independently by the author in [2] and Wheelon in [3] and that allows 
one to find exact values for a large class of convergent series of rational terms. The 
method has been extended to many series in the papers of Lesko and Smith [3j and 
Efthimiou [5]. 

We will outline both the author's and Wheelon's variations of the method in 
order to compare the similarities and differences, and then use them to compute the 
series ([1]). Finally, we apply the methods to additional series that generalize ([1]). 
The results thus obtained seem to be new. 

The Core Idea of the Method Let a series of the form Ylkei '^k^k^ where / is 
a subset of Z, be given. Then it is convenient to write only one of the factors, say 
Vk, as an integral transform 



Problem 819. For n > 1, evaluate 




(1) 




and thus 
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Assuming that the order of the operations of summation and integration can be 
exchanged 



k&I \k&I 

In this article we shall always exchange the order of the two operations assuming 
that the reader knows how to reason for its validity. Details on this may be found 
in the original papers [2] and [3] where the Laplace transform is used. If one can 
find an explicit function h{t) = Y2kei'^kVk{t), then he has succeeded in writing the 
initial series in a simple integral representation: 

b 



^UkVk = / h{t)f{t)dt. 
kei 



If, furthermore, the integration can be performed, then analytic answers for the 
initial series are obtained. 



Series of rational terms Consider a series for which / = {1,2,...}, nfc = l, and 

Vk = Q{k)/P{k), where Q{k) and P{k) are two polynomials in k: 

To ensure convergence, we may assume that degQ + 2 <degP. 

Problem 819 of The College Mathematics Journal is clearly a special case of ^ 
lor Q{k) = I and P{k) = kik + l){k + 2) ■ ■ ■ {k + n). 



The author's variation In the method presented in [2], the expansion of the 
general rational term of the series to partial fractions plays a central role. We shall 
then assume that the sum S can be written in terms of partial fractions: 

oo e rrii . 



ii=ij=i ^ ^ 

where the constants Aij are uniquely determined by the partial fraction decompo- 
sition of each summand. Details for this decomposition can be found in [2]. 



Using the identity 



1 



(L - 1)! 



-oo 

^L-lg-Ax^^^ (3) 
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we write the series in integral form (which is vahd only if aj > —1, for all i): 

S = EEEtAv / x-^e-C^-^)^^^. (4) 

fc = l 2 = 1 j = l 

By interchanging the integration over x with the addition over k, we find 

i=l j=2 ^ 

+ y Ail dx . 

i=i 

To further simplify this sum we need to know the values of the constants Aij. The 
integrals may be expressed in terms of the polygamma functions [2j. 



Andreoli's problem In problem 819 of The College Mathematics Journal the 
rational function 1/P{k) can be easily expanded in partial fractions 

^ A- 



with 

(-1)* fn 



This can be proved easily by multiplying equation ^ hy k + j and then taking the 
limit k ^ —j. We thus have 



°o 1 1 e"^ ^fn 



= — / —{l-e-'')'^dx 

n\ Jq 1 - e ^ 

1 f + co 

= — / e-^(l-e-^)"^i(ix . (6) 
nl Jo 

By a change of variables u = e~^, 

(l-u)"-idn = . (7) 



°° I 1 fi 

Ep(^ = 



n n\ 
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A generalization One can easily generalize Andreoli's proposal: For n > 1, we 
can evaluate 



"^^"'^^ £ [a + kb][a + {k + l)b][a + {k + 2)b]---[a + {k + n)b] 
Following the method described above, we find 

1 r+oo 

The change of variables u = gives the integral the form 

Sia,b) = ^j\'^^'-\l-ur-'du. 
The last integral is a B-function integral: 



(8) 



{l-uy-^u'^-^du = jB{j,n) , m,n,l>^. (9) 



The sum S{a, b) is thus equal to B{a/b + 1, n) which value can be easily computed 
through the F-function. Finally, we get 



1 " 1 



nb a + ib 

i=l 



Wheelon's variation The method presented by Wheelon [3] for the computation 
of the series ([2]) assumes that Q{k) = 1 and P{k) can be factorized as a product 
of of first order monomials. In this method, the convolution theorem as used by 
Feynman plays a central role. In particular, starting from 

1 dx 



ab Jo [ax + 6(1 — x)]'^ ' 
one can show the general formula 

-L i=i 1=1 i=i 



i=l i=l 



dx\dx2 ■ ■ ■ dxi- 



1 • 



In physics literature, this formula is known as Feynman's integral and the x's are 
called Feynman's parameters. 



5 



If Q{k) = 1 and P{k) given by 

P{k) = {k + ai)™^ {k + a2T^ ...{k + atT^ , (10) 

where the Oj, i = 1,2,... ,i are distinct real numbers, none of them a negative 
integer and all rm are positive integers, then the series ([2]) is written 

r{Em) ^ (n|^r"')(i-|:^.r-' 

S = -j—^ ^ / /•■■/ — 1 dx\dx2- ■ ■ dx(,-x 



-1 £-1 



n r(mi) k=i u 



i=l < E < 1 

i = l 



1 



^ e-1 
y e 

„ i-i i-i 

■■■ dxidX2---dxl.l(Y[xT'~'){l-^Xir^-' X 

n r(mi) " ti 
i=i < E < 1 

i = l 

, < 1 

e-1 

e -y[ai+l+J2 ^ii'^i-ae)] 

-oo (J]m,-1) e «=i 

^ 1-e-y 

To pass from the first expression to the second, we used formula 1^. Wheelon does 
not present explicitly the above formula; he only describes his method in a simple 
case and then he lists exact results for series for which P{k) is a product of up to 
four distinct monomials. The result for Andreoli's problem (which will be computed 
in the next section) does not appear in [3j. 



Andreoli's problem again When nii = 1, for all i, Feynman's integral can be 
simplified considerably. In particular, one can write 

(^ — 1)! / dxi / dx2--- / dxe^i X 



aia2 ■■■ai 



\a\ + xi(a2 - ai) H h X£-2(a£-i - (^1-2) + 2;£_i(a<> - a<?_i)]^ 
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If £ = n + 1 and = k + i — 1, i = 1, 2, . . . , n + 1, then 
1 



k{k + l){k + 2)---{k + n) 



and 



^ A;(A; + l)(A; + 2)- 



1 /-Xi rl^n-l 

n! / dxi / (ix2 ■ • • / dxn 

JO JO 



[A: + E 



/ dxi / dx2-- - I dXn :^ 

Jo Jo Jo [fc+^^.]n+l 



i=l 



1 



— / dxi dx2-- - I dXn I y"e 
Jo Jo Jo Jo 



1 f+°° „ -y[k+j:x,] 

dxr, I V e dy 



n\ JO 



e-y^{l-e~yTdy 







^_^k{k + l){k + 2) ■ ■ ■ {k + n) 



1 r+oo 
1 r+oo 

{l-e-yr~^e~y dy , 



nl Jo 

which is identical to equation ([6]) found by the partial fractions method. 



(11) 



A Related Result Equation (jlip is valid for any positive number k (we shall 
rename it x), not necessarily an integer. If instead we sum over n 

^x(x + l)(x + 2)---(x + n) Jq ^ 



n=0 



el u"" du , 



where we have made a change of variable u = e ^ . In this expression we can now 
expand e~" in a Taylor series and perform the elementary integrals to find 



oo -. 

V - 

^ x(x + + 2) • • • (x + n) 

n=0 



111 11 

+ ^ 



X l!x + l 2!x + 2 



(13) 



Additional Results The series proposed by Andreoli and its generalization ([8]) 
are well-known: see page 234 of [7]. Equations (fT3]) and (fT2]) are also found in [7] 
(page 317). Equation ([13]) — although not equation (fT2]) — is presented as an exercise 
(ex. 126, page 272) in Knopp's book f8]. In this section we extend the results 
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presented above to other series. The results obtained seem to be new; we have not 
found them in [7] or [8] . 

Using the Feynman integral and equation ([3]), we can show that 



1 



1 



^ / u^-^ (1 - u'y-' du . (14) 



k{k + 2){k + A)---{k + 2n) ni 
By renaming k as x, summing over n, and following the steps that led to equation 
(fT3l). we can find 



n=0 



x{x + 2){x + 4) • ■ • (x + 2n) 



/ u--^e-^'/^du 



,1/2 



111 11 

+ 



x 21! x + 2 2^2! x + 4 
On the other had, summing in (jl4p over k gives: 



. (15) 



1 



1 



(1 -n^)" du . 

1 — u 



^ k{k + 2){k + 4)---{k + 2n) 2"n! Jq 

We could have arrived at the same integral representation if we had used the partial 
fraction decomposition. The integral of the right hand side may be rewritten as 



(l-n")"-i(l + u)dn = / {l-u')''-'du+ {l-uT''udu. 

Jo Jo 

These integrals in turn can be computecj^ by the use of the B-function ([9]). Finally, 

we find 



y I = i , 

^^k{k + 2){k + A)---{k + 2n) 2nV(2n-l)!! (2n) 



1 1 

+ 



(16) 



Empowered with the results that have been presented so far, one can ask if we 
can prove similar results for the series 



n=0 

oo 

sii) = 



1 



x{x + i){x + 2£) ■ ■ ■ {x + ni) 
1 

^J^ k{k + £){k + 2i) ■ ■ ■ {k + nl) 



Indeed, it is straightforward to show that 

-1 

F{x 







111 11 

+ 



X il\ x + i im x + 2i 



(17) 



^Actually, the second integral is basic. 
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and that 



S{£) 




1-u 



1 



du . 



The integral of the right hand side may be rewritten as 



^0 



/ (l-n^)"-i(l + n + u2 + ... + -u' 



) du 



that is, a sum of integrals of the form Q. We thus find 



Sii) = ^( 
in \ 



1 1 1 



+ ••• + 



{ln-{l-\)Y 



1 




where the symbol !^ is an extension of the double factorial and stands for the product 
of all numbers between 1 and In — m that, when divided by give residue i — m. For 



example, (3n)!!! = (Sn)!^ = 3-6-9 • • • (3n), (3n-l)!!! = (Sn-l)!^ = 2-5-8 • • • (3n-l), 



Conclusion Taking a problem proposed by Andreoli as the origin of this article, 
we have attempted to draw the attention of the reader to the use of integral trans- 
forms in the computation of exact values for series and make him appreciate the 
ease and transparency of the method. In the process, we have derived some results 
([T71) and (fT8|) that seem to be new or, at least, do not appear in the commonly used 
references [7] and [5]. 

The method can be applied to a wide variety of series as explained in the articles 
of Efthimiou [21 [5] and Lesko and Smith [3j and the book of Wheelon [3] . The reader 
is encouraged to look at these references for additional information. Hopefully, he 
will use the method to search for new results. He may begin with the series: 



Acknowledgements Few days before the posting of the present paper on the 
archives, I received a message from M. L. Glasser drawing my attention to one 
of his papers [9] in which he had extended Wheelon's method as described in his 
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with [2], [3] and [5] although the presentation differs from these papers. 



and so on. 
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